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Abstract 


The  objective  of  this  investigation  is  to  determine  to  what  extent  a 
highly-maneuverable  aircraft's  overall  performance  capability  is 
enhanced  by  thrust-vectoring  nozzles.  The  resulting  performance 
capabilities  are  compared  to  a  baseline  configuration  with  non-thrust- 
vectorlng  nozzles  to  determine  the  effects  and  advantages  of  thrust 
vectoring. 

The  results  indicate  that  the  use  of  vectored  thrust  can 
significantly  increase  an  aircraft's  performance  capability  in  turning 
flight.  The  greater  the  demand  on  the  aircraft  in  a  turning  combat 
scenario,  the  more  the  aircraft  utilizes  its  thrust-vectoring  capability 
to  complete  the  task.  The  results  also  indicate  that  the  use  of 
vectored  thrust  in  other  phases  of  flight  —  such  as  cruise, 
acceleration  and  climb  —  only  slightly  increases  an  aircraft's 
performance  capability. 
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I .  Introduction 


Background 

A  fighter  aircraft's  performance  capability  mainly  consists  of  a 
combat  mission  radius  associated  with  its  close-in  combat 
maneuverability.  By  evaluating  the  aircraft's  performance 
characteristics  in  segments  which,  when  added  together,  define  the 
combat  mission  profile,  the  design  of  the  fighter  aircraft  can  be  broken 
down  into  three  parts:  climb,  cruise  and  combat.  The  aircraft's 
maximum  climb  performance  is  determined  by  the  engine's  maximum 
available  thrust;  optimum  cruise  performance  is  determined  by  its 
aerodynamic  and  propulsive  efficiency;  and  combat  performance  is 
determined  by  maximum  turning  capabilities  —  dictated  by  structural  and 
angle  of  attack  limits. 

In  an  effort  to  improve  a  fighter  aircraft's  performance  without 
changing  the  basic  design,  the  Air  Force  is  currently  investigating  the 
use  of  vectored  engine  thrust  by  developing  the  F-15  short  takeoff  and 
landing  (STOL)  demonstrator.  By  installing  advanced  aerodynamic  and 
propulsion  control  concepts  on  an  "off-the-shelf"  F-15,  the  Air  Force 
expects  the  F-15  STOL/maneuver  technology  aircraft  to  "yield  combat 
performance  and  runway  flexibility  improvements  of  a  magnitude  normally 
associated  with  the  development  of  an  entirely  new  aircraft”  (1).  The 
primary  modifications  to  the  F-15B  two-seat  flight  test  aircraft  will  be 
"the  addition  of  controllable  canards  above  the  engine  inlets  forward  of 
the  main  wing  of  the  aircraft”  and  "two-dimensional  nozzles  at  the  rear 
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of  the  F-15's  Pratt  and  Whitney  F100  engines.”  The  F-15  STOL 
demonstrator  is  expected  to  begin  flying  in  1988,  testing  these 
"potential  applications  to  the  next  generation  of  Air  Force  combat 
aircraft."  The  Advanced  Tactical  Fighter,  nicknamed  as  a 
"superfighter,"  might  employ  vectored  engine  thrust  —  "an  especially 
welcome  feature  in  a  superfighter  because  of  its  small  wings"  (2). 

Several  studies  present  optimal  controls  and  trajectories  to 
minimize  an  aircraft's  turning  time,  one  aspect  of  an  aircraft's  overall 
performance  capability.  Humphreys,  Hennig,  Bolding  and  Helgeson  (3) 
investigated  three-dimensional  aircraft  dynamics.  Johnson  (4)  included 
the  possibility  of  in-flight  thrust  reversal.  Finnerty  (5)  constrained 
the  maneuver  to  the  vertical  plane.  Brinson  (6)  considered  the  effects 
of  sideforce  and  Schneider  (7)  considered  the  effects  of  vectored  thrust 
in  reducing  the  time  to  turn.  Johnson,  Brinson  and  Schneider  all  found 
substantial  benefits  in  utilizing  additional  controls  to  reduce  an 
aircraft's  turning  time  in  three-dimensional  flight.  This  study,  in 
contrast,  investigates  a  highly-maneuverable  tactical  aircraft's  overall 
performance  capability  enhanced  by  thrust-vectoring  nozzles  and 
restricts  turning  flight  to  the  horizontal  plane. 

Problem  Statement 

This  investigation  seeks  the  thrust  vector  angle  schedules  which 
will  optimize  various  air  combat  energy  maneuverability  parameters  — 
including  specific  excess  power,  maximum  sustained  and  Instantaneous 
turn  rate,  and  the  power  setting  required  to  sustain  a  given  load  factor 


and  energy  level. 


The  controls  to  be  optimized  are  nozzle  deflection  angle  and  thrust 
sideslip  angle  (where  applicable).  These  controls  are  to  be  considered 


unconstrained  initially  and  subsequently  constrained  to  practical 
physical  limits  as  dictated  by  the  F-15  STOL  demonstrator  design. 

To  evaluate  the  effects  of  thrust  vectoring  on  fighter  aircraft 
performance,  the  optimal  thrust  vector  angle  schedules  and  associated 
performance  results  will  be  obtained  and  compared  against  the  results  of 
a  non-thrust-vectoring  configuration. 

Assumptions 

Several  assumptions  regarding  the  aircraft,  its  dynamics  and  the 
controls  are  incorporated  into  this  study.  These  assumptions  are  not 
only  necessary  to  reduce  the  complexity  of  the  problem  to  manageable 
proportions,  but  are  common  in  this  type  of  study  as  a  source  of 
meaningful  comparison  of  results.  The  assumptions  are: 

1.  the  aircraft  is  a  point  mass 

2.  flat,  non-rotating  earth 

3.  NASA  1962  Standard  Atmosphere  (8) 

4.  constant  gravitational  acceleration 

5.  coordinated  turn  (no  sideforce) 

6 •  instantaneous  controls 

Summary  of  Current  Knowledge 

Aircraft  performance  analyses  done  within  the  Air  Force  and  industry 
have  traditionally  been  done  using  graphical  and  analytical  methods  with 
underlying  assumptions  that  the  angle  of  attack  and  thrust  vector  angle 
are  small  and  consequently  neglected.  In  the  case  of  highly- 


maneuverable  aircraft  utilizing  large  angles  of  attack  and  high  thrust- 
to-weight  ratios  (such  as  the  F-15  and  F-16),  the  negligible  angle  of 
attack  assumption  becomes  less  accurate  but  will  yield  acceptable 
conservative  results.  In  the  case  of  the  F-15  STOL  demonstrator  where 
thrust-vectoring  nozzles  will  be  used,  the  thrust  vector  angle  becomes  a 
primary  control  variable  and  will  directly  influence  the  canard 
deflection. 

Current  aircraft  performance  computer  programs  in  use  throughout  the 
Air  Force  are  typically  generic  —  meaning  that  they  are  applicable  to 
nearly  every  type  of  Air  Force  aircraft,  such  as  cargo,  transport, 
trainer,  attack  and  fighter  aircraft.  To  keep  the  programs  from 
becoming  too  complex,  the  small  angle  of  attack  and  thrust  vector  angle 
assumptions  were  built  into  the  logic  since  there  was  no  real  need  for 
these  capabilities  at  the  time.  With  the  advent  of  the  F-15  STOL 
demonstrator  program  and  its  eventual  Advanced  Tactical  Fighter 
applications,  the  inclusion  of  angle  of  attack  and  thrust  vector  angle 
as  additional  control  variables  in  optimal  performance  analyses  has 
become  a  necessity  for  aircraft  performance  engineers. 

There  is  only  one  aircraft  currently  in  production  and  service  in 
the  United  States  that  utilizes  thrust-vectoring  nozzles  —  the  AV-8 
Harrier,  flown  by  the  United  States  Marine  Corps.  The  Harrier  is 
capable  of  vectoring  its  four  exhaust  nozzles  through  angles  of 
approximately  100  degrees,  giving  it  vertical  and  short  takeoff  and 
landing  capability  (V/STOL).  Low  speed  pitch,  roll  and  yaw  attitudes 
are  controlled  by  six  small  reaction  control  jets  placed  in  appropriate 
locations  on  the  aircraft  and  are  integrated  into  the  flight  control 


system.  The  Harrier  uses  its  thrust-vectoring  capabilities  mainly  to 
provide  direct  lift  since  the  exhaust  nozzles  are  symmetrically  located 
about  the  aircraft's  center  of  gravity.  The  F-15  STOL  demonstrator,  on 
the  other  hand,  uses  its  thrust-vectoring  capabilities  to  provide  both 
direct  lift  and  pitch  control  through  its  Integrated  control  system  that 
governs  the  nozzle  deflection  angle/canard  deflection  interaction. 

Approach 

Definitions  of  all  aspects  of  the  aircraft  performance  problem 
including  models  and  constraints  are  discussed  in  Section  II.  The 
derivations  of  the  equations  of  motion  of  an  aircraft  in  straight  and 
level  flight,  climbing  flight  and  turning  flight  are  presented  in 
Section  III.  These  equations  of  motion  include  the  angle  of  attack  and 
thrust  vector  angle  contributions  which  are  traditionally  neglected.  A 
direct  comparison  of  the  traditional  equations  of  motion  is  also  covered 
in  this  section.  The  changes  necessary  to  transform  the  traditional 
equations  of  motion  in  an  existing  aircraft  performance  computer  program 
to  those  that  include  angle  of  attack  and  thrust  vector  angle  are 
discussed  in  Section  IV.  The  methods  used  in  optimizing  the  aircraft's 
performance  with  vectored  thrust  are  presented  in  Section  V.  The 
results  of  this  study  are  then  discussed  and  the  use  of  thrust  vectoring 


is  compared  to  the  non-thrust-vectoring  baseline  in  Section  VI. 
Finally,  conclusions  and  recommendations  are  given  in  Section  VII. 


II.  The  Aircraft  Performance  Problem 


Before  results  can  be  analyzed  and  compared,  It  Is  necessary  to 
completely  define  all  aspects  of  the  problem.  The  problem  will  be 
defined  In  terms  of  the  different  phases  of  flight  to  be  Investigated, 
the  characteristics  which  model  the  aircraft,  atmospheric  properties, 
and  practical  physical  constraints  on  the  control  variables. 

Phases  of  Flight 

The  different  phases  of  flight  are  determined  by  setting  certain 
state  variables  and/or  their  derivatives  equal  to  zero  and  letting  the 
other  state  variables  vary  within  the  optimization.  The  phases  of 
flight  Included  In  this  study  are: 

1 .  straight  and  level  flight 

a.  acceleration 

b.  cruise 

2.  climbing  flight 

3.  turning  flight 

a.  sustained  turn 

b.  instantaneous  turn 

The  takeoff  and  landing  phases  are  not  included  since  they  require  an 
optimal  trajectory  analysis  beyond  the  scope  of  this  study.  This  study 
concentrates  on  the  aircraft's  mission  performance  applications,  which 
treat  the  takeoff  and  landing  portions  of  a  mission  profile  as  segments 
of  flight  that  use  fuel  with  no  time  or  range  considerations.  With  the 


phases  of  flight  Included  in  this  study  and  appropriate  approximations 
for  takeoff  and  landing,  a  typical  fighter  aircraft  mission  profile  can 
be  analytically  performed. 

Equations  of  Motion 

The  equations  of  motion  for  flight  of  a  point  mass  aircraft  over  a 
flat,  non-rotating  earth  are  derived  by  Miele  (9:42-49)  as 


X  =•  V  cos  Y  cos  X  (1) 

Y  =  V  cosY  sinx  (2) 

h  »  V  sinY  (3) 

mV  *  T  cos£  cos^  -  D  -  mg  sinY  (4) 

X  cosM  cosY  -  YsinM  »  ^  {T  cose  sinv  -  Q  +  mg  sinp  cosy)  (5) 

'  '  1 

X  sinu  cosY  +  YcosU  *  {T  sin£  +  L  -  mg  cosP  cosY^  (6) 


The  definitions  of  all  the  variables  used  in  this  study  are  contained  in 
the  List  of  Symbols  preceding  the  Abstract  and  are  shown  pictorially  in 
Figure  1.  Since  this  study  does  not  allow  for  sideforce,  Q  *  0. 
Rearranging,  the  equations  of  motion  become 


X  -  V  cosy  cos x 


(7) 


Y  ■  V  cos  Y  sin  X 


( 


h  ■  V  sin  Y 


"  T  D 

V  *  g  (y  cos  e  cos  v  -  —  -  sinY  ^ 


V  cos  y 


T  L 

{^(cose  sinv  cosy  +  sine  siny)  +  ^  siny) 


o  x  L 

Y  =  v  {y  (sinE  cosy  -  cose  sin^  siny)  +  ^  cosy  -  cosy) 


These  equations  are  written  in  the  wind  axes  and  describe  the  aircraft 
motion  with  respect  to  an  earth-fixed  coordinate  frame.  The  state 
variables  are  X,  Y,  h,  V,  X  and  Y.  The  variables  y,  e  and  v  are 
control  variables.  The  aircraft  weight  (W)  is  considered  constant  for 
the  point  performance  calculations,  but  will  vary  when  used  in  quasi¬ 
steady  integrations. 

The  lift,  drag  and  thrust  forces  will  be  discussed  in  the  next 
paragraphs  and  will  give  rise  to  three  more  control  variables. 


Aircraft  Characteristics 


The  common  coefficient  forms  of  the  aerodynamic  lift  and  drag  forces 

are 


L  - 


%>V  S  CL 


(13) 


D  =  ^PV  S  C  D 


(14) 


From  incompressible  aerodynamic  and  thin  airfoil  theories,  the  lift  and 
drag  coefficients  can  be  expressed  as 


°L 


CT  +  c,  =  CT  +  (L  0.  +  r  +  r  6 

TJB  c  0tra  a  L0  L6 

WB  WB  c  c 


(15) 


WB  c 


(16) 


which  shows  that  the  lift  coefficient  is  a  linear  function  of  angle  of 
attack  (a)  and  canard  deflection  angle  (£c)  and  that  the  drag 
coefficient  varies  parabolically  with  lift  coefficient.  These  simple 
yet  descriptive  relationships  give  rise  to  the  terms  "linear  lift  curve 
slope”  and  "parabolic  drag  polar"  and  identify  the  control  variables 
associated  with  the  aircraft's  aerodynamics.  The  aerodynamic  lift  and 
drag  data  used  to  model  the  aircraft  in  this  study,  however,  have  been 
derived  from  flight  test  and  analytical  methods  and  exhibit  the  same 
trends  as  the  linear  lift  curve  slope  and  parabolic  drag  polar.  These 


data,  contained  in  Appendix  A,  realistically  model  the  aircraft's 
aerodynamics  with  angle  of  attack  and  canard  deflection  angle  as  control 
variables . 

The  engine  data  is  also  derived  from  flight  test  and  analytical 
methods  and  does  not  resemble  the  ideal  jet  theory  where  engine  thrust 
is  considered  only  as  a  function  of  altitude.  The  thrust  and  fuel  flow 
characteristics  can  be  written  ideally  as 

T-T»ax*  <17>  afc-iUffl  +  tQ)  (18) 


where  Tm  is  a  function  of  velocity  and  altitude  and  the  thrust  power 
setting  (tt)  is  the  control  variable  associated  with  the  aircraft's 
propulsion  system.  These  data  are  also  contained  in  Appendix  A. 

The  formulation  of  the  lift,  drag  and  thrust  relationships  has 
introduced  several  aircraft  parameters  and  three  control  variables:  a  ,  6 
and  x  .  Approximate  values  for  the  aircraft  parameters  in  the  case  of 
the  F-15  STOL  demonstrator  at  a  Mach  number  of  0.8  and  an  altitude  of 
30,000  feet  are: 
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S  =  608  ft 


W  =  35,000  lb 


j 


Tmax  =  10.900  lbT  (Military  Power) 


i 


T _  -  22,968  lbT 


max 


(Maximum  Afterburner) 


t^£  =  0.833  lbf/(lb,p  hr)  (Military  Power  and  below) 


t^£  =  2.142  lbj/(lb.j,  hr)  (Maximum  Afterburner) 


tQ  =  1000  lbf/hr 


These  values  were  used  initially  with  the  linear  lift  curve  slope, 
parabolic  drag  polar  and  engine  thrust  equations  to  test  the  techniques 
developed  for  optimizing  aircraft  performance  with  thrust  vector  angle, 
while  the  results  presented  in  this  study  are  based  on  the  data  in 
Appendix  A. 

Atmospheric  Model 

The  NASA  1962  Standard  Atmosphere  (8)  was  used  for  this  study  and  is 
limited  to  the  first  two  layers  of  the  lower  atmosphere.  This 
atmosphere  is  modelled  from  sea  level  to  the  tropopause  (36,089  feet)  as 
a  linearly  decreasing  temperature  layer,  and  above  the  tropopause  as  an 
isothermal  layer  —  the  troposphere  and  stratosphere,  respectively.  The 
aircraft  performance  program  used  in  this  study  is  programmed  with  the 
complete  set  of  equations  to  model  the  NASA  1962  Standard  Atmosphere. 


Control  Variable  Constraints 

Three  control  variables  —  the  throttle  setting  (tt),  angle  of  attack 
(ot)  and  canard  deflection  (<5  )  —  are  constrained  by  physical 
considerations.  The  thrust  can  neither  be  greater  than  the  maximum 
thrust  nor  less  than  the  minimum  thrust  of  the  particular  power  setting 
being  used.  The  F-15  is  capable  of  throttling  both  the  military  and 
afterburner  power  settings.  The  minimum  thrust  for  the  military  power 
setting  is  taken  to  be  zero,  while  the  minimum  thrust  for  the 
afterburner  power  setting  is  taken  to  be  100%  military  power.  The  scope 
of  this  study,  however,  will  only  be  concerned  with  throttling  military 
and  not  the  afterburner  power  setting  for  cruise  and  turning  performance 
calculations.  Since  the  throttle  setting  is  defined  in  equation  (17)  in 
terms  of  the  maximum  thrust,  the  throttle  setting  is  limited  to 

0  1  for  military  power 

^  =  1  for  afterburner  power 

Each  lifting  surface  is  limited  by  its  own  aerodynamic  lift  limit 
and  the  total  lift  generated  is  limited  by  the  aircraft's  maximum  load 
factor.  Individually,  the  angle  of  attack  and  canard  deflection  can  not 
exceed  their  stall  points;  collectively,  the  angle  of  attack  and  canard 
deflection  can  not  generate  more  lift  than  is  tolerated  by  the 
aircraft's  maximum  load  factor.  Since  the  canard  is  used  to  trim  out 
the  pitching  moment  caused  by  the  deflection  of  the  thrust  vector,  the 
canard  surface  is  assumed  to  be  designed  such  that  its  stall  point  is 
never  reached  in  normal  hrust  vectoring  applications  to  maintain  pitch 
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control  authority.  In  this  case,  the  aircraft  is  constrained  by  an 

angle  of  attack  limit  (defined  by  a  CL  of  1.6)  and  a  thrust  vector 

max 

angle  limit  and  not  constrained  by  a  canard  deflection  limit. 

In  turning  performance,  the  velocity  where  the  aircraft  is  at  its 
maximum  angle  of  attack  and  maximum  load  factor  is  the  corner  velocity. 
The  corner  velocity  is  the  velocity  at  which  the  aircraft  achieves  its 
maximum  turn  rate.  At  speeds  below  the  corner  velocity,  the  aircraft  is 
limited  by  its  maximum  angle  of  attack.  At  speeds  above  the  corner 
velocity,  the  aircraft  is  limited  by  its  maximum  load  factor  (9  g's). 

No  constraints  were  placed  on  the  thrust  angle  of  attack  (E)  and 
thrust  sideslip  angle  (v).  While  the  F-15  STOL  demonstrator  will  have 
limits  on  its  nozzle  deflection  angle  (^) 

0  jp  o 

-20  1  5n1+20 

this  angle  was  allowed  full  range  in  order  to  determine  how  much  range 
of  thrust  vectoring  would  be  exploited  if  it  were  available.  The 
effects  of  constraining  the  two  thrust  vector  angles  were  later  examined 
and  are  discussed  in  Section  VII. 

The  aircraft's  velocity  is  constrained  by  the  lesser  of  the 

following:  a  dynamic  pressure  limit  determined  by  the  aircraft's 

2 

structural  and  flutter  characteristics  (q  **  2131  lb/ft  );  or  the 
aircraft's  Mach  number  limit  determined  by  temperature  or  engine  limits 
(M  ■  2.5).  No  constraint  was  placed  on  the  bank  angle  (p). 


Ill .  Aircraft  Performance  Equations  of  Motion 

Miele's  equations  of  motion  for  flight  of  a  point  mass  aircraft  over 
a  flat,  non-rotating  earth  as  discussed  in  Section  II  use  the  thrust 
sideslip  angle  (v)  and  thrust  angle  of  attack  (c)  as  the  angles  of 
successive  rotation  to  which  the  wind  axes  must  be  subjected  in  order  to 
turn  the  velocity  vector  in  a  direction  parallel  to  the  thrust  vector. 
When  a  pair  of  body-fixed  reference  lines  are  introduced  into  the 
derivation  of  the  equations  of  motion  (as  depicted  in  Figure  2),  the 
thrust  angle  of  attack  control  variable  is  found  to  be  dependent  upon 
the  aircraft's  angle  of  attack  (a)  —  a  control  variable  derived  from 
the  aircraft's  aerodynamic  characteristics. 


Figure  2  —  Definition  of  Thrust  Vector  Angle 
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The  nozzle  deflection  angle  (5N)  replaces  the  thrust  angle  of  attack  as 
the  second  angle  of  rotation  In  the  definition  of  the  thrust  vector 
orientation,  with  the  body-fixed  engine  centerline  replacing  the 
velocity  vector  as  the  reference  orientation.  Since  this  study  allows 
for  no  sideslip,  the  thrust  sideslip  angle  definition  is  unchanged.  The 
nozzle  deflection  angle  and  thrust  angle  of  attack  are  related  by  the 
equality 

£  *  5N  +  T  +  “  -  (*9) 

The  equations  of  motion  now  become 


X  -  V  cosy  cos  x  (20) 

Y  =  V  cos  y  sin  x  (21) 

h  -  V  sinY  (22) 

V  -  g  {  —  cos(SN  +  t  +  a)  cos  v  -  ^  -  sinY^  (23) 


■  a  X 

X  55  V  cos  Y  {wtco8(5N  +  T  +  a)  sinV  cos  ^  +  sin(5N  +  T  +  a)  siny  ] 

+  ~  sin  y  }  (24) 

Y  -  ft  J[sin(6  +  x  +  a)  cosy  -  cos(6  +  t  +  a)  sinv  siny] 


The  first  three  equations  are  kinematic  relationships  and  the  second 
three  equations  are  dynamic  relationships. 


Trimmed  Flight 

Since  the  aircraft  is  considered  to  be  a  point  mass  with  all  forces 
acting  at  the  center  of  gravity,  there  must  be  no  net  moments  acting  on 
the  aircraft.  In  the  case  of  the  F-15  STOL  demonstrator,  however,  a 
deflection  of  the  nozzle  will  cause  a  moment  since  the  nozzle's  point  of 
action  is  not  at  the  center  of  gravity.  The  addition  of  the  canard  in 
the  F-15  STOL  demonstrator  design  allows  for  balancing  the  moment 
induced  by  the  nozzle  deflection  by  deflecting  the  canard  appropriately. 
This  interaction  results  in  an  additional  constraint  equation  relating 
the  canard  deflection  and  the  nozzle  deflection.  Since  the  F-15  STOL 
demonstrator  design  does  not  have  a  sideforce  generator  to  balance  a 
moment  caused  by  deflecting  the  nozzle  out  of  the  aircraft's  plane  of 
symmetry,  the  thrust  sideslip  angle  (v)  must  be  set  to  zero  for  trimmed 
flight. 

The  additional  constraint  equation  relates  the  canard's  aerodynamic 
characteristics  to  the  thrust  (T),  nozzle  deflection  angle,  the  aircraft 

angle  of  attack  and  the  aircraft's  dimensional  characteristics  (1  , 

x 

lc  ,  1 N  ,  lty)  by  balancing  the  moments  induced  by  the  nozzle  and  canard 

Z  X  z 

deflections : 


lc  [Lc  cos®  +  Dc  8ina]  +  lc  [Lc  sin<*  -  Dc  cos**] 

x  z 

+  [T  sin(  <5m  +t  )]  +  1M  [T  cos(Sn+t)]  -  0 


deflection  control  variable  in  the  same  manner  as  the  lift  and  drag  of 
the  aircraft  are  functions  of  the  angle  of  attack  control  variable. 


Functional  Relationships 

The  final  set  of  equations  governing  the  motion  of  the  point  mass 
aircraft  is  composed  of  three  kinematic  relationships,  three  dynamic 
relationships  and  one  trim  constraint  equation.  The  trim  constraint 
equation  transforms  the  aerodynamic  characteristics  of  the  aircraft  from 
a  function  of  four  variables  to  a  function  of  six  variables: 

L  ■  L(h,  V,  a  ,  6  )  trim  L  *  L(h,  V,  a  ,  <5  ,  6  ,  ir  ) 

c  _  C  N 

«=> 

D  -  D(h,  V,  a  ,  6  )  constraint  D  *  D(h,  V,  a  ,  Sc  ,  ,  it  ) 

The  kinematic  and  dynamic  equations  have  the  following  functional 


relationships : 


* ) 


X  =■  X(V,  Y,  x  ) 

Y  -  Y(V,  Y,  X) 

h  =  h(V ,  y) 

Y  =  V(h,  V,  Y,  a,6c  ,&N  , 

x  =  x(h,  V,  Y,  a,6c  ,6  ,  it,  y) 

Y  -  Y(h,  V,  Y,  a,6c  , 6^  ,  v ,  y) 

The  aircraft  performance  computer  program  used  in  this  study 
utilizes  point  performance  optimization,  treating  the  kinematic 
equations  as  integrals  of  motion  and  the  dynamic  equations  as  the 
differential  system.  To  accomplish  this,  the  aircraft's  downrange  (X) 
and  crossrange  (Y)  locations  along  with  its  altitude  (h)  and  velocity 
(V)  must  be  known  prior  to  optimization.  Since  the  program  already  has 
routines  to  optimize  the  various  aircraft  performance  parameters  with 
altitude  and  velocity,  the  optimization  of  the  aircraft's  performance 
with  thrust  vectoring  can  be  considered  as  a  point  sub-optimization 
problem,  or  an  optimization  within  an  optimization. 

Of  the  twelve  state  and  control  variables,  seven  remain  unspecified 
and  are  considered  dependent  variables.  Since  there  are  four  equations 
in  the  differential  system  (including  the  trim  constraint  equation), 
there  are  three  degrees  of  freedom  in  which  to  optimize  the  aircraft's 
performance  (9:  48-51).  By  restricting  flight  to  the  vertical  and 
horizontal  planes  and  defining  a  phase  of  flight  by  constraining  certain 
variables  within  that  plane,  the  number  of  degrees  of  freedom  in  which 
to  optimize  the  aircraft's  performance  can  be  reduced  significantly. 
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Flight  In  the  Vertical  Plane 

This  category  of  trajectories  is  represented  by  a  vertical  plane  of 
symmetry  which  implies 

Y  »  constant  (27 

V  =  0  (28 

Equation  (27)  implies  that  Y  -  0,  and  equation  (21)  implies  that 
X  *  X  =  0.  The  equations  of  motion  become 


X  =  V  cos  Y  (29 

h  *  V  sin  Y  (30 

V  =  J  [T  cos(6N  +  t+oi)-D-w  sinY]  (31 

Y  “  _g_  [T  8in( 6  +x+a)+L-W  cosy]  (32 

WV 

X  >  0  (33 

0  ■  1  [Lc  cosa  +  Dc  sina]+  1  [L  sina  -  D  cosa] 

x  z 

+  1N  [T  sin(6N+  x)]  +  1N  [T  cos(5N +x  )]  (34 

X  z 


Rewriting,  the  remaining  equations  of  motion  are 


X  =  V  cos  Y 


h  =  V  sin  Y 


V  [T  cos(SN  +t+a)-D-W  siny] 


y  -  [T  sin(  5N+T+a)+L-W  cosy] 


0  =  lc  [Lc  cosa  +  Dc  sinct]  +  [Lc  sina  -  Dc  cosa] 
x  z 

+  1N  [T  sin(6N+T  )]  +1N  [T  cos(6n  +  t)1 

X  Z  ’’ 


This  reduced  system  of  equations  for  flight  in  the  vertical  plane  has 
two  kinematic  relationships  and  three  dynamic  equations.  Of  the 
seven  dependent  variables  previously  unspecified  for  point  performance 
optimization,  five  variables  —  "»  6n»  a  .  6C  and  Y  —  remain 

unspecified,  yielding  two  degrees  of  freedom. 

One  common  application  of  flight  in  the  vertical  plane  is  that  of 
straight  and  level  flight,  where  the  flight  path  angle  (Y)  and  its 
derivative  (Y)  are  set  to  zero.  This  implies  that  the  aircraft  will 
maintain  a  constant  altitude  during  this  flight  phase.  The  equations  of 
motion  reduce  to 


X  -  V 


V  -  J  [T  cos(  SN+  t+o  )  - 


>V, 


0  »  T  sin(<5M  +  t+ci)+l-W 


0*1  [L  cosa  +  D  sina]  +1  [L  sina  -  D  cosa] 

c  c  c  *  c  c  c 

x  z 

+  1N  [T  sin(6N  +  T)]  +  1N  [T  cos(5n  +t  )] 

X  z 


with  dependent  variables w  ,  6^  ,  a  and  6^.  Two  cases  of  straight  and 
level  flight  are  of  primary  interest  to  the  aircraft  performance 
engineer  —  acceleration  and  cruise.  Although  the  equations  of  motion 
are  the  same,  the  techniques  involved  in  finding  the  aircraft's  optimal 
acceleration  and  cruise  controls  are  entirely  different. 

In  the  case  of  maximizing  an  aircraft's  acceleration,  the  throttle 
setting  control  variable  is  set  at  some  constant  value  (normally  at  the 
maximum  —  tt»  l  for  maximum  acceleration).  This  reduces  the  number  of 
dependent  variables  to  three  (  6N,  a  and<5c  ),  but  since  the  performance 
parameter  to  be  optimized  (V)  is  solved  for  explicitly  in  the  first 
dynamic  equation  (41),  the  number  of  equations  in  the  differential 
system  is  also  reduced  by  one.  The  resultant  differential  system  has 
three  dependent  variables  and  two  equations,  leaving  just  one  degree  of 
freedom  in  optimizing  the  aircraft's  acceleration  capability.  Since 
this  study  is  interested  in  finding  the  effects  of  thrust  vectoring  on 
the  performance  of  an  aircraft,  the  one  control  variable  chosen  to  be 
the  degree  of  freedom  is  the  nozzle  deflection  angle  ($N). 

In  the  constant  altitude  cruise  case,  the  throttle  setting  (tt)  is 
the  parameter  that  is  to  be  minimized.  Since  the  cruise  condition  is  a 
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non-accelerating  flight  phase,  the  aircraft  s  acceleration  (V)  must  be 
set  to  zero.  The  throttle  setting  appears  explicitly  in  all  three 
dynamic  equations  but  also  appears  implicitly  in  the  first  two  as  a 
functional  dependency.  Thus  the  differential  system  has  four  dependent 
variables  and  three  equations,  leaving  just  one  degree  of  freedom. 

Again  the  nozzle  deflection  angle  is  the  control  variable  that  is  chosen 
to  be  the  single  degree  of  freedom. 

The  traditional  equations  of  motion  for  level  flight  in  the  vertical 
plane  neglect  the  angle  of  attack  and  thrust  vector  angle  contributions: 


X  -  V 


V  -  §  [T  -  D] 


L  -  W 


Since  the  aircraft  weight  is  known,  the  lift,  drag  and  subsequently  the 
aircraft's  acceleration  or  cruise  throttle  setting  are  easily  calculated 
—  results  from  a  system  of  equations  with  no  degrees  of  freedom. 

Another  common  application  of  flight  in  the  vertical  plane  is  that 
of  climbing  flight,  where  the  parameter  to  be  optimized  is  the  time  rate 
of  change  of  the  altitude  or  climb  rate  (h).  The  climb  rate  is  solved 
for  explicitly  in  the  second  kinematic  equation  but  can  be  inserted  into 
the  first  dynamic  equation  by  a  direct  substitution.  The  equations  of 
motion  reduce  to 


.-•V-  .'-V  V  .'-'.V. .  v  V-V  V.V-'A'i 
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I 


X  *  V  cos  y 


(47) 


h  *  77  [T  cos(6  +t+a)-D--V] 
W  N  g 


(48) 


Y  *  [T  sin(6N  +r+a)+L-W  cosy] 


(49) 


0*1  [L  cosa  +  D  sinal  +1  [L  sina  -  D  cosal 

c  c  c  c  c  c 

X  z 


+  1N  [T  8ln(6N  +T  )]  +  1  [T  cos(6n  +t  )] 

X  L  z 


(50) 


An  alternate  way  of  expressing  the  climb  rate  Is 


h  * 


[T  cos  (6N  +  t  +  a)  -  D]  v 


1  +  V  dV 
8  dh 


W 


(51) 


V  dV 


where  —  —  Is  the  acceleration  correction  term.  For  climbing  flight  at 


a  constant  true  airspeed,  the  acceleration  terms  are  zero: 


or 


V  w  * 

V  ■  V 

dV  ' 

wfl  - 

8  V  *g 

dh  h 

w  ‘ 

V 

dV 

-v*o 

and  g 

dh 

(52) 


(53) 


The  equations  of  motion  for  non-accelerating  climbing  flight  become 


X  *  V  cos  y 


(54) 
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(55) 


h  =  |  [T  cos(6n  +  x  +  a)  -  Dj 

Y  =  WV  [T  sin^<SN  +  T  +  a)+L-w  cosy]  (56) 

0  =  1  [L  cos  a+  D  sinal  +1  [L  sina  -  D  cosal 

C  C  C  J  c  c  c 

X  z 

+  1N  [T  sin(6N+T)]  +  1N  [T  cos(6N  +T  )]  (57) 

X  z 

The  climb  rate  equation  (55)  is  Identical  to  the  equation  for  the 
aircraft's  level  acceleration  (41)  when  multiplied  by  a  constant  factor. 
Since  the  non-accelerating  climb  rate  and  level  acceleration  are  the 
parameters  to  be  optimized  in  their  respective  cases,  their 
optimizations  should  result  in  the  same  values  for  the  control 
variables . 

This  relationship  between  the  constant  altitude  acceleration  and  the 
non-accelerating  climb  is  the  basis  for  the  specific  energy  concept 
described  by  Rutowski  (10).  In  his  approach  to  the  performance  problem 
using  energy  concepts,  the  author  computes  the  time  rate  of  change  of 
the  aircraft's  energy  in  terms  of  its  climb  rate  and  acceleration: 

^(§)  =Ps  =  h+  |v=|[T  cos(6N  +  t  +  a)  -  D]  (58) 

V  ' 

For  the  constant  altitude  case,  Pg  =  —  V.  For  the  non-accelerating 
case,  Pg  *  h.  Optimizing  the  specific  excess  power  (Pg)  optimizes  the 
accelerating  climb  rate  or  the  varying  altitude  acceleration  providing 
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the  desired  simplification  with  no  loss  in  accuracy. 


Since  the  optimization  of  the  specific  excess  power  can  be  done  by 
the  techniques  developed  for  the  level  acceleration,  the  dynamic 
equation  for  the  time  rate  of  change  of  the  flight  path  angle  (y)  is  not 
used  when  optimizing  climb  rate.  This  parameter  can  be  calculated, 
however,  for  the  non-accelerating  climb  by  making  use  of  the 
relationship  between  the  climb  rate  and  the  flight  path  angle: 


h  ■  V  sinY 


(59) 


I 

I 

I 

I 

Flight  in  the  Horizontal  Plane 

This  category  of  trajectories  is  represented  by  the  mathematical 
condition 

i 

i 

I 

h  =  constant  (63) 

I 

I 

I 

I 

Equation  (63)  implies  that  h  =  0,  and  equation  (30)  implies  that 

Y  =  y  m  o.  The  remaining  equations  of  motion  are  I 


X  =  V  cos  X  (64) 

Y  =  V  sin  X  (65) 

V  =  g  [T  cos(6N  +  t  +  a)  -  D]  (66) 

X  =  ^  [T  sin(6N  +  t  +  a)  sinP  +  L  sinM]  (67) 

0  =  T  sin(<5^  +  t  +  a)  cosM  +  L  cosM  -  W  (68) 


0  =  lc  [Lc  cos'*  +  Dc  sina]+  lc  [Lc  sin<*  -  Dc  cos'*] 
x  z 

+  1N  [T  sln(5N  +  x)]  +  lN  [T  cos(6n  +  t)]  (69) 

X  z 

This  reduced  system  of  equations  for  flight  in  the  horizontal  plane  has 
two  kinematic  relationships  and  four  dynamic  equations.  Of  the  seven 


dependent  variables  previously  unspecified  for  point  performance 


optimization,  six  variables  —  ,  <5C  ,  u  and  X  —  remain 

unspecified.  The  heading  angle  (X)  is  only  used  in  the  kinematic 
relationships  and  is  not  included  in  the  dynamic  equations,  so  the 
differential  system  has  only  five  dependent  variables,  yielding  one 
degree  of  freedom. 

A  common  application  of  flight  in  the  horizontal  plane  is  that  of 
the  instantaneous  turn,  where  the  aircraft  is  accelerating  or 
decelerating  (V  y  0).  This  implies  that  the  aircraft  will  maintain  a 
prescribed  turn  rate  at  the  specified  throttle  setting  but  will  not  turn 
at  a  constant  altitude,  constant  velocity  combination  (energy  state) 
during  this  flight  phase.  Since  this  study  restricts  the  aircraft  to 
turn  only  in  the  horizontal  plane,  the  altitude  must  remain  constant  and 
any  change  from  the  initial  energy  state  will  be  reflected  in  a  change 
in  velocity.  By  maximizing  the  aircraft's  acceleration  (or  minimizing 
its  deceleration),  the  aircraft's  energy  gain  (or  loss)  is  optimized. 

The  performance  parameter  to  be  optimized  (V)  is  solved  for  explicitly 
in  the  first  dynamic  equation  (66),  so  the  number  of  equations  in  the 
differential  system  is  reduced  by  one.  Since  the  throttle  setting  is 
specified,  this  leaves  one  degree  of  freedom:  the  nozzle  deflection 
angle  (  6^. 

Another  application  of  flight  in  the  horizontal  plane  is  that  of  the 
sustained  turn,  where  the  aircraft  is  not  accelerating  or  decelerating 
(V  =  0).  This  implies  that  the  aircraft  will  maintain  a  constant 
altitude  and  velocity  during  the  flight  phase.  The  equations  of  motion 


reduce  to 


X  =*  V  cos  X 


Y  =  V  sinX 


=  T  cos(  fiN+  t  +  <*  )  -  D 


[T  sin(  <5^  +  t  +a  )  sin y  +  l  sinP ] 


0  =>  T  sin(  6N  +  T  +  a)  cos  y  +  L  cos  U  -  W 


0  =  1  [L  cos  a  +  D  sina]  +1  [L  sina  -  D  cos«] 


c  1  c 
x 


c  1  c 
z 


•>  1N  [T  8in(6N+T  )]  +  1N  [T  cos(5N  )] 


with  dependent  variables  n,  f a  ,  <$c  and  u  .  There  are  two  cases  of  the 
sustained  turn  flight  phase  that  are  of  primary  interest  to  the  aircraft 
performance  engineer  —  the  maximum  sustained  turn  rate  and  the  minimum 
throttle  setting  required  to  maintain  a  specified  sustained  turn  rate. 
Although  the  equations  of  motion  are  the  same,  the  techniques  involved 
in  finding  the  maximum  sustained  turn  rate  and  the  minimum  throttle 
setting  are  entirely  different. 

In  the  case  of  maximizing  the  aircraft's  sustained  turn  rate,  the 
throttle  setting  control  variable  is  set  at  some  constant  value 
(normally  at  the  maximum  —  it  *  1).  This  reduces  the  number  of 
dependent  variables  to  four  a,  <5C  and  y),  but  since  the  performance 
parameter  to  be  optimized  (x)  is  solved  for  explicitly  in  the  second 
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dynamic  equation  (73),  the  number  of  equations  in  the  differential 

I 

system  is  also  reduced  by  one.  The  resultant  differential  system  has 
one  degree  of  freedom  and  again  will  be  chosen  to  be  the  nozzle 
deflection  angle. 

In  the  case  of  finding  the  minimum  throttle  setting  required  to  j 

t 

maintain  a  specified  sustained  turn  rate,  the  throttle  setting  appears 

t 

t 

explicitly  in  all  four  dynamic  equations  but  also  appears  implicitly  in 

! 

the  first  three  as  a  functional  dependency.  Thus  the  differential 
system  has  five  dependent  variables  and  four  equations,  and  will  use  the 
same  control  variable  as  the  degree  of  freedom  as  the  other  turn  cases. 

The  traditional  equations  for  turning  flight  in  the  horizontal  plane 
neglect  the  angle  of  attack  and  thrust  vector  angle  contributions: 

n-|  (76) 

X  =  V  cos  X  (77) 


Y  -  V  sin  X 


(78) 


V 


w  [T  -  D] 


(79) 


X 


_S_ 

wv 


L  sinM 


(80) 


For  a  sustained  turn  the  thrust  is  equal  to  the  drag,  which  determines 
the  load  factor  through  the  drag  polar  and  in  turn  determines  the  turn 
rate.  For  an  instantaneous  turn  the  specified  turn  rate  determines  the 
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load  factor,  which  la  turn  determines  the  drag  and  the  aircraft's 


acceleration. 


Summary 


Restricting  flight  to  the  vertical  and  horizontal  planes  greatly 
reduces  the  magnitude  and  complexity  of  the  aircraft  performance 
problem.  Applications  of  this  study  reduces  the  flight  In  the 
vertical  and  horizontal  planes  to  optimization  problems  with  one  degree 
of  freedom. 
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IV.  Aircraft  Performance  Computer  Program 

NSEG  (11),  a  segmented  aircraft  performance  analysis  program  widely 
in  use  throughout  the  Air  Force,  uses  the  traditional  equations  of 
motion  to  calculate  an  aircraft’s  performance  capability.  Implementing 
the  full  equations  of  motion  discussed  in  Section  III  would  require  a 
major  overhaul  of  the  existing  equations  which  would  necessitate  a 
complete  checkout  for  all  applications.  On  the  other  hand,  adjusting 
the  thrust  and  drag  data  to  account  for  the  changes  in  the  equations  of 
motion  would  leave  the  existing  equations  intact  while  only  changing  the 
data  lookup  routines.  Currently,  the  aircraft's  aerodynamic  drag  data 
is  input  as  a  function  of  the  required  lift  force,  Mach  number, 
altitude,  drag  index  and  center  of  gravity.  The  engine’s  propulsion 
characteristics  (net  thrust  and  fuel  flow  data)  are  input  as  a  function 
of  power  setting,  Mach  number  and  altitude. 

The  NSEG  program  utilizes  a  thrust-drag  accounting  system  that 
treats  the  airframe  and  throttle-dependent  drags  (ram,  spillage  and 
nozzle  drag)  as  separate  and  unrelated.  This  system,  when  coupled  with 
the  small  angle  of  attack  and  thrust  vector  angle  assumptions,  reduces 
Miele's  equations  of  motion  (20-25)  to  the  traditional  equations  of 
motion  (60-62  and  76-80),  where  the  thrust  (T)  is  actually  the  net 
thrust  (F^)  defined  by 

F  =■  F  -  D  -  D  ...  -D  .  (81) 

N  g  ram  spillage  nozzle 

and  is  the  resultant  propulsive  effort  of  the  engine.  The  thrust  and 
drag  appear  in  the  traditional  equations  of  motion  only  as  the  quantity 
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(T  -  D),  which  renders  the  thrust-drag  accounting  system  arbitrary  as 
long  as  it  remains  consistent. 

This  thrust-drag  accounting  system  becomes  very  important,  however, 
when  the  complete  equations  of  motion  are  used.  In  this  case,  the 
thrust  (T)  is  the  engine's  gross  thrust  (F  ),  which  is  the  actual  force 

o 

being  exerted  through  the  nozzle.  The  total  drag  (D)  will  then  have  to 
include  the  throttle-dependent  drags  as  discussed  above.  As  the  gross 
thrust  is  vectored  through  a  positive  angle  (SN+  T  +  <*),  there  is  a 
loss  in  thrust  corresponding  to  a  gain  in  lift.  To  illustrate  this 
point,  part  of  Figure  1  is  repeated  here  with  additional  definitions. 


Figure  4  —  Thrust's  Contribution  to  Lift 

When  (  6n  +  t  +a)  ■  0,  the  engine's  propulsive  force  is  just  Fg.  When 
(<5^  +  t  +  a)  >  0,  the  propulsive  force  along  the  flight  path  is 
Fg  cosC  +  x  +  cx ),  The  loss  in  thrust  along  the  flight  path  is  then 
the  difference  between  these  two  forces: 


(82) 


v‘- 

v* 

I 


_y 
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F  -  Fg  [1  -  cos(5n  +  t  +a  )] 

In  flight,  this  loss  in  thrust  may  be  expressed  as  an  increase  in  drag. 
In  drag  coefficient  form,  this  increase  is  given  by 

F  [1  -  cos  (6  +  t  +  ot)] 

AC  -  -£ - ^ -  (83) 

D  h  p  v2  s 
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The  corresponding  increase  in  lift  due  to  vectoring  the  thrust,  as  shown 
in  Figure  4,  is 

&L  -  F  sin(<5  +  x  +a  )  (84) 

g  N 


s 

.  1 

.v 
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In  lift  coefficient  form,  this  increase  is  given  by 


F  sin  (6  +  t  +  a) 

4Cl  ■  J - r - 

p  V  s 


(85) 


The  total  lift  and  drag  can  then  be  expressed  as 


+  AC, 


aero 


(86) 


C 


D 


CD  +  ACD 
aero 


(87) 


where  is  the  lift  generated  by  the  aircraft's  aerodynamics  and  CQ 

aero  aero 

is  its  corresponding  drag.  For  v  -  0,  NSEG's  traditional  equations  of 

motion  and  thrust-drag  accounting  system  defined  by  equation  (81),  along 

with  these  adjustments  to  the  aerodynamic  data  will  provide  the  same 
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results  as  If  the  complete  set  of  equations  of  motion  were  reprogrammed 
into  NSEG.  The  required  data  now  includes  the  aircraft's  aerodynamic 
data  and  the  engine's  net  thrust  and  fuel  flow  characteristics  as 
before,  as  well  as  the  aircraft's  lift  curve  slopes  and  the  engine's 
gross  thrust  characteristics.  The  program  uses  the  gross  thrust  data  to 
adjust  the  aeodynamic  data  accordingly  and  uses  the  net  thrust  as 
before.  For  v  4  0,  the  turn  rate  equation  (24)  and  the  time  rate  of 
change  of  the  flight  path  equation  (25)  each  have  an  additional  term  to 
be  considered  and  programmed.  Since  v  has  been  determined  to  be  zero 
for  this  aircraft  in  Section  III,  the  v-  0  case  has  not  been  addressed. 

Included  in  the  programming  changes  to  account  for  thrust  vectoring 
is  the  trim  constraint  equation  (26).  The  moment  induced  by  vectoring 
the  thrust  is  counteracted  by  inducing  an  equal  and  opposite  moment  with 
the  canard.  With  the  canard  deflection,  additional  lift  and  drag  are 

added  and  taken  into  account  in  C  L  and  CD  respectively.  The 

aero  aero 

thrust  (T)  in  the  trim  equation  is  actually  the  gross  thrust  (Fg )  being 
exerted  through  the  nozzle  that  creates  the  moment  about  the  center  of 


gravity. 


V.  Optimization  Techniques 


As  discussed  in  Section  III,  the  aircraft  performance  problem  can  be 
broken  down  into  several  distinct  phases  of  flight,  each  of  which  are 
one  degree  of  freedom  optimization  problems.  Since  the  functional  form 
of  the  performance  parameters  to  be  optimized  are  not  explicitly  known, 
numerical  techniques  must  be  used  to  search  for  the  optimal  control 
values.  Two  techniques  —  equal  interval  search  and  quadratic 
interpolation  —  were  coupled  together  to  provide  an  efficient  and 
effective  optimization  algorithm  with  an  acceptable  interval  of 
uncertainty. 

The  equal  Interval  search  technique  was  used  to  find  the  interval  in 
which  the  optimal  solution  is  located.  In  this  procedure,  the  control 
variable  chosen  to  be  the  degree  of  freedom  is  set  at  some  lower  bound 
and  is  incremented  until  an  upper  bound  is  met  or  until  an  optimal 
Interval  has  been  detected.  At  each  point,  the  performance  parameter  to 
be  optimized  is  calculated  and  compared  to  its  value  at  the  previous 
point.  The  optimal  interval  is  detected  when  the  performance  parameter 
calculated  is  not  the  optimal  value  when  compared  to  its  previous  value. 
The  optimal  Interval  is  then  defined  as  the  interval  bounded  by  the 
points  on  either  side  of  the  optimal  solution  determined  by  the  equal 
interval  search. 

The  quadratic  interpolation  technique  uses  the  three  points  that 
define  the  optimal  interval  along  with  their  functional  values  to 
determine  the  coefficients  of  a  quadratic  equation  that  approximates  the 
performance  parameter  function  in  this  interval.  By  elementary 
calculus,  the  optimal  control  value  can  then  be  determined  analytically 


and  the  optimal  performance  parameter  value  can  be  calculated  directly. 
Using  a  second  order  Lagrangian  interpolation  polynomial,  the  quadratic 
equation  can  be  represented  by 

F(x)  ■  ax^  +  bx  +  c  **  F. (x)  +  F9(x)  +  F9(x)  (88) 


where 


x-x2 

X~X3 

xrx2 

xrx3 

x-x1 

X-X3 

X2“X1 

X2'X3 

X-Xj^ 

x-x2 

X3’X1 

x3‘x2 

f  (xL) 


f(x2) 


Fl(x) 


F9(x)  - 


F3(x) 


Equating  coefficients,  the  quadratic  equation’s  coefficients  are  found 
to  be 


a  -  (x3-x2^  f^xp  +  ^xi-x3^  f^x2^  +  (X2“XP  f(x3) 

(X3-X1)  (x  -x  )  (x  -x  ) 


(x2+x3) (x2~x3)f (Xj)  +  (x1+x3)(x3-x1)f(x2)  +  (x1+x2) (x1~x2>f (x3> 
(x^j-x^)  (x3~x2)  (x^Xj) 


The  quadratic  equation  is  then  differentiated  once  and  set  to  zero,  and 
the  optimal  control  value  is  easily  determined: 


2 

F(x)  =  ax  +  bx  +  c 


(95) 


F  (x)  =  2ax  +  b  =  0 


(96) 


Equation  (89)  implies  that 


x 


_b_ 

opt  “  2a 


(97) 


Even  though  this  optimization  algorithm  is  not  the  best  algorithm 
that  could  be  used  to  find  a  function's  optimal  value,  it  was  used 
because  the  NSEG  program  already  utilizes  the  quadratic  interpolation 
technique  in  optimizing  simple  functions.  In  these  existing 
applications,  NSEG  systematically  picks  three  points  in  a  given 
interval,  calculates  their  functional  values,  and  "curvefits"  the 
function  to  find  an  approximate  optimal  solution.  The  equal  interval 
search  and  quadratic  interpolation  techniques  were  coupled  together  in 
this  new  application  so  that  a  fairly  accurate  optimal  solution  would 


result . 


VI.  Results 


Since  the  objective  of  this  investigation  is  to  determine  to  what 
extent  a  highly-maneuverable  aircraft's  overall  performance  capability 
is  enhanced  by  thrust-vectoring  nozzles,  the  results  are  shown  against  a 
baseline  configuration  with  non-thrust-vectoring  nozzles  such  that  a 
direct  comparison  can  be  made.  The  results  are  shown  for  three  cases: 

Case  1:  Results  from  using  the  traditional  equations 
of  motion  with  angle  of  attack  and  thrust 
vector  angle  neglected 

Case  2:  Results  from  using  the  complete  equations  of 
motion  with  no  thrust  vectoring 

Case  3:  Results  from  using  the  complete  equations  of 

motion  at  the  optimum  nozzle  deflection  angle. 

With  the  thrust-vectoring  nozzles  being  used  in  flight  to  enhance 
the  aircraft's  combat  maneuverability,  the  results  are  shown  for  various 
points  throughout  the  air  battle  arena.  Parrott  (12)  defines  the  air 
battle  arena  as  the  part  of  the  fighter  aircraft's  envelope  in  which 
most  of  its  close-in  combat  maneuvering  takes  place.  The  air  battle 
arena  definition  assumes  a  threat  engagement  at  a  high  velocity,  high 
altitude  condition.  As  the  dogfight  develops,  the  aircraft  dissipates 
its  energy  by  maneuvering  outside  of  its  sustained  envelope  and 
eventually  ends  up  at  a  low  velocity,  low  altitude  condition.  In  order 
to  be  effective  in  a  dogfight,  the  aircraft  must  be  able  to  maneuver 
without  excess  energy  loss  and  must  be  able  to  disengage  from  the 


opponent  at  any  time.  This  requires  high  sustained  load  factors  and 
high  level  flight  acceleration  capabilities  throughout  the  air  battle 
arena. 

Figures  5  and  6  illustrate  the  F-15  STOL's  air  battle  arena 
superimposed  on  the  aircraft's  sustained  envelope  for  various  load 
factors.  Tables  1  through  9  show  the  results  for  the  critical  Mach- 
altitude  points  of  the  air  battle  arena  for  the  cases  previously 
discussed.  Appendix  8  contains  graphical  data  for  the  first  case  and 
Appendix  C  contains  graphical  data  for  the  second  and  third  cases. 
Appendix  D  contains  graphical  data  for  subsonic,  high  altitude,  long 
range  cruise  conditions  for  mission  range  comparisons.  The  units  for 
the  different  aircraft  performance  parameters  are  as  follows: 


Specific  Excess  Power  (Ps)  :  ft/sec 

Sustained  Load  Factor  :  g's 

o 

Sustained  Turn  Rate  :  /sec 

Thrust  Required  :  lb 

Specific  Range  :  NM/lbf 


Figures  27  through  34  show  that  an  aircraft's  sustained  load  factor 
and  turn  rate  are  significantly  improved  through  the  use  of  thrust 
vectoring,  especially  at  low  altitudes  and  low  velocities.  The  tables 
of  results  quantitatively  show  an  increase  in  sustained  load  factor  as 

much  as  0.3  g  at  maximum  A/B  power  and  0.23  g  at  military  power,  and  an 

o 

Increase  in  sustained  turn  rate  as  much  as  1.3  /sec  at  maximum  A/B  and 

o 

0.5  /sec  at  military  power  in  the  air  battle  arena.  Figures  35  through 


42  show  that  an  aircraft's  specific  excess  power  is  increased 
substantially  through  the  use  of  thrust  vectoring,  especially  at  higher 
load  factors.  The  aircraft's  specific  excess  power  is  shown  unchanged 
at  level  1  g  flight.  Figures  43  through  46  show  that  lower  thrust  is 
required  to  sustain  a  specified  load  factor  with  thrust  vectoring, 
especially  at  lower  speeds  that  are  not  part  of  the  air  battle  arena. 
Figures  50  through  52  show  that  the  aircraft's  high  altitude,  subsonic 
cruise  performance  is  improved  slightly  with  thrust  vectoring. 
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TABLE  I 

Results : 

M  =  0.7/10,000 

feet 

Case  1 

Case 

2  Case  3 

Military  Power 

Ps  6  Ig 

266 

267 

267  <?  0° 

Ps  S  3g 

188 

190 

194  @  6° 

Ps  «  5g 

-48 

-30 

9  @  14° 

Sustained  Load  Factor 

4.74 

4.83 

5.04  @  14° 

Sustained  Turn  Rate 

11.3 

11.6 

12.1  @  14° 

%  Thrust  Required  for 

3g  Sustained  Turn 

55.0 

54.1 

53.5  @  1° 

Maximum  A/B  Power 

ps  @  ig 

687 

687 

687  @  0° 

ps  0  5g 

372 

401 

442  @  8° 

Ps  0  9g 

-1255 

-1064 

-678  @  25° 

Sustained  Load  Factor 

6.28 

6.47 

7.10  @  21° 

Sustained  Turn  Rate 

15.2 

15.6 

17.2  @  21° 
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TABLE  II 


Results:  M  =  0.8/10,000  feet 


Case  1 

Case  2 

Case  3 

Military  Power 

Ps  «  is 

259 

259 

259  @  0° 

ps  0  3g 

201 

202 

204  @  3° 

Ps  0  5g 

51 

56 

65  <3  8° 

Sustained  Load  Factor 

5.38 

5.45 

5.65  @  10' 

Sustained  Turn  Rate 

11.3 

11.5 

11.9  (3  10' 

%  Thrust  Required  for 

57.6 

57.5 

57.1  @  4° 

3g  Sustained  Turn 


Maximum  A/B  Power 

ps  0  ig 

794 

794 

794  @  0° 

Ps  6  5g 

587 

594 

606  @  5° 

Ps  0  9g 

-595 

-484 

-169  @  21 

Sustained  Load  Factor 

7.59 

7.78 

8.43  <3  19' 

Sustained  Turn  Rate 

16.1 

16.5 

17.9  <3  19* 
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TABLE  III 


Results:  M  =  0.9/10,000  feet 


Case  1 

Case  2 

Case  3 

Military  Power 

Ps  0  lg 

232 

232 

232 

0  0° 

Ps  0  3g 

181 

183 

184 

0  2° 

Ps  @  5g 

58 

61 

67 

0  e 

Sustained  Load  Factor 

5.55 

5.61 

5.75 

0  9° 

Sustained  Turn  Rate 

10.4 

10.5 

10.8 

0  9° 

%  Thrust  Required  for 

3g  Sustained  Turn 

67.0 

66.6 

66.4 

@  f 

Maximum  A/B  Power 

Ps  @  lg 

906 

906 

906 

0  0° 

Ps  0  5g 

732 

737 

747 

0  5° 

Ps  @  9g 

-40 

19 

164 

0  15° 

Sustained  Load  Factor 

8.86 

9.00 

9.00 

0  17° 

17.0  0  17° 


Sustained  Turn  Rate 


16.7 


17.0 


TABLE  IV 


Results:  M  =  1.0/10,000  feet 


Case  1 

Case  2 

Case  3 

Military  Power 

Ps  @  ^ 

-6 

-6 

-5  0  1° 

Ps  @  3g 

-96 

-94 

-86  0  6° 

Ps  @  5g 

-283 

-278 

-252  @  10° 

Sustained  Load  Factor 

- 

- 

- 

Sustained  Turn  Rate 

- 

- 

- 

%  Thrust  Required  for 

3g  Sustained  Turn 

— 

— 

Maximum  A/B  Power 

Ps  @  lg 

837 

837 

838  0  1° 

Ps  @  5g 

560 

569 

603  0  7° 

Ps  @  9g 

-291 

-248 

-65  @  15° 

Sustained  Load  Factor 

7.98 

8.13 

8.73  0  16° 

14.9  0  16° 


Sustained  Turn  Rate 


13.6 


13.8 


TABLE  V 

Results : 

M  =  0.8/20,000 

feet 

Case  1 

Case  2 

Case  3 

Military  Power 

Ps  @  lg 

216 

216 

216  0  0° 

Ps  @  3g 

117 

119 

123  0  7° 

Ps  0  5g 

-266 

-239 

-158  0  23 

Sustained  Load  Factor 

3.93 

3.99 

4.17  0  12 

Sustained  Turn  Rate 

8.5 

8.6 

9.0  0  12 

X  Thrust  Required  for 

67.8 

67.3 

66.4  0  7° 

3g  Sustained  Turn 

Maximum  A/B  Power 

Ps  @  lg 

603 

604 

604  @  1° 

Ps  @  5g 

121 

163 

260  0  14 

Ps  @  9g 

-2473 

-2068 

-1403  0  29 

Sustained  Load  Factor 

5.34 

5.50 

6.00  @  18 

Sustained  Turn  Rate 

11.7 

12.0 

13.2  0  18 

TABLE  VI 


Results:  M  =  1.0/20,000  feet 


* 

Case  1 

Case  2 

Case  3 

Military  Power 

•  \ 

p, «  u 

51 

51 

52  e  1° 

„  1 

Ps  9  3g 

-80 

-77 

-67  @  9° 

l.'. 

Ps  «  5g 

-399 

-386 

-337  @  14° 

8 

Jfl 

Sustained  Load  Factor 

2.03 

2.05 

2.12  @  7° 

A' 

Sustained  Turn  Rate 

3.1 

3.2 

3.3  6  7° 

■r.' 

■r: 

%  Thrust  Required  for 

3g  Sustained  Turn 

- 

- 

- 

4 

i 

Maximum  A/B  Power 

Ps  @  ^ 

664 

665 

665  0  0° 

1 

Ps  @  5g 

214 

236 

300  @  12° 

Pg  @  9g 

-1099 

-995 

-789  0  21° 

>4 

V 

Sustained  Load  Factor 

5.82 

5.95 

6.44  @  17° 

Sustained  Turn  Rate 

10.2 

10.4 

11.3  0  17° 

58 


JVuV\.V 


Wff 


mmmm 


TABLE  VII 


Results:  M  =  1.2/20,000  feet 


Case  1 

Case  2 

Case  3 

Military  Power 

pse  i8 

-525 

-525 

-525  0  0° 

Ps  9  3* 

-539 

-537 

-534  @  5° 

Ps8  5g 

-815 

-806 

-763  0  13 

Sustained  Load  Factor 

- 

- 

- 

Sustained  Turn  Rate 

- 

- 

- 

X  Thrust  Required  for 

3g  Sustained  Turn 

- 

- 

- 

Maximum  A/B  Power 

P  @  lg 
s  6 

410 

410 

410  0  0° 

P  @  5g 
s  G 

121 

136 

192  0  8° 

Ps  @  9g 

-1262 

-1175 

-889  0  21* 

Sustained  Load  Factor 

5.50 

5.63 

6.07  0  12 4 

Sustained  Turn  Rate 


8.0 


8.2 


8.9  0  12 


TABLE  VIII 


Results:  M  =  0.9/30,000  feet 


Case  1 

Case  2 

Case  3 

Military  Power 

P  @  lg 

s 

144 

144 

144 

@  o° 

P  @  3g 

s 

-18 

-12 

7 

@  10° 

P  @  5g 

s 

-698 

-651 

-474 

0  27° 

Sustained  Load  Factor 

2.89 

2.94 

3.05 

e  io° 

Sustained  Turn  Rate 

5.6 

5.7 

5.9 

0  10° 

%  Thrust  Required  for 

- 

- 

99.9 

0  io° 

3g  Sustained  Turn 


Maximum  A/B  Power 


Ps  @  lg 

479 

480 

480 

0 

0° 

PS  @  5g 

-362 

-286 

-92 

0 

20 

Ps  @  9g 

-5592 

-4777 

-3092 

0 

48 

Sustained 

Load  Factor 

4.23 

4.34 

4.71 

0 

19 

Sustained 

Turn  Rate 

8.5 

8.7 

9.5 

0 

19 

TABLE  IX 


Results:  M  =  1.4/30,000  feet 


Military  Power 


Case  1 


A  -  .  ■r.  -  .  -  . 


Case  2 


Case  3 


rs  @  lg 

-592 

-592 

-592  @  0° 

Ps  @  3g 

-701 

-698 

-691  @  7° 

/ 

Ps  @  5g 

-1197 

-1177 

-1078  @  14* 

Sustained  Load  Factor 

- 

- 

- 

Sustained  Turn  Rate 

- 

- 

- 

%  Thrust  Required  for 

— 

_ 

_ 

.A 

3g  Sustained  Turn 

*•  - 

i 

Maximum  A/B  Power 

Ps  @  lg 

378 

378 

378  @  0  c 

Ps  @  5g 

-228 

-191 

-74  @  13 

(1 

Ps  @  9g 

-2366 

-2215 

-1733  @  23 

> 

Sustained  Load  Factor 

4.27 

4.36 

4.71  @  12 

‘  j 

Sustained  Turn  Rate 

5.5 

5.6 

6.1  @  12 

v.v. 


v'/.'A-', 


VII.  Conclusions  and  Recommendations 


Two  major  conclusions  are  drawn  based  upon  the  results  of  this  study: 

1.  Thrust-vectoring  nozzles  can  substantially  increase  a  fighter 
aircraft's  turning  performance  throughout  the  air  battle  arena. 
The  +  20°  nozzle  deflection  angle  constraint  would  limit 

the  aircraft's  optimal  performance  in  only  a  few  instances. 

2.  Level  1  g  flight  does  not  benefit  from  the  use  of  vectored  thrust 
since  the  aircraft's  wing  loading  is  small  compared  to  that  of 
turning  flight. 

Two  recommendations  are  offered  for  future  research: 

1.  This  study  concentrated  on  the  air-to-air  combat  applications  of 
a  fighter  aircraft's  performance.  Air-to-ground  applications 
should  be  investigated  for  the  ordnance  delivery  scenario. 

2.  Investigate  the  class  of  aircraft  that  might  utilize  thrust 
sideslip  vectoring  and  sideforce  generation  and  the  type  of 
trajectories  that  may  result. 
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REYNOLD'S  NUMBER  CORRECTIONS 


I  NO AO  1 1 
I  A01X«< 


l  A01Y-: 
ATABOl  = 


•.oooa 

.000/ 


>. .  10000 


i  »  L  •  0  «  L  . 
-.0012  ,- 


•.0009*  - 
-.0008  ,- 


-.0008,- 
-.0007  ,- 


-.0007,- 
■  • 000  o »  — 


-.0006,- 

'.0006 »  — 


■ . 0006  ,  - 
■•0005  »  — 


. .20000.  , 


2. L. 5. 1.6 

• G0C6. • 0 . 


.0005. .0. 
. 00  04. . 0 . 


•0004.. 0. 
.0004. «0. 


.0004. . 0  . 

. 00  0  3. .0. 


.0003. .0. 
. 00 03. . 0 . 


. 000  3. .0  . 
. 0002 ..0. 


300 00.  .40000.  t  50000. .60000.  .70000, 


fl.8.2.0f2.2fc.5f 

.0009, .0020, .0033, .0047,  ' 


.0007. . 0015,. 0025, . C 
.0006  ,  .  00  14  ,.  002  3  ,.  <; 


.0006, . Ou 13, .0022, .9 031,  S  M  *  1  •  0 
.0005, .0012,. 0  02  Of. 0029,  SM-1.2 


.0005 , .0012,. 00 l 9,. 00 2 
.0005,  .0011, .Odd, .002 


.0004, .0010, .0017, .0023,  Irt-l.d 
.0004,  .0010,. 0016,. 0022,  tfi-2.0 


.CG03..0009..0 
.0003, . 0008  .  . 0 


=  2.2 

=  2.5 


F-15  FLIGHT  TEST  DRAG  POLARS 


INDA02 


I  A02X- 
I A02Y- 


AT  AB02 


.0228, . 
.0493. 


.2105, 

.  r'  2  2  2  ,  . 


50 

63 


.0220, . 

. 0453 , 


.2203, 
.0217, . 


.0437, 
. 2290, 


.0215,. 

.0445, 


.2 

. 0216, . 


. 0466  , 
.2400, 


=  1  . 


27, 
21 , 


*.0,  . 

7  5  ,  •  8 


.  4  ,  . 
1.3, 


0225, 

.0580 


.2400 

0215, 


.  05  30 
.2465 


0213, 

.0523 


.  2525 
0213. 


.0520 
.  2643 


0213, 

.0534 


.2700 

0213, 


.05 

.27 


.9 


6  ,  •  7  ,  .  3 
1.4, 1.5 


.0230,  . 
,.0685, 


,.2715, 

.0220,. 


,.0635, 
,.  2  79  3, 


.02  15,. 
,.0625, 


, .2880, 
.02 16,  . 


, .0633, 
,.3000, 


.0220,. 

,.0646, 


15,  .2 
,  .95  , 


,  .  8  5  , 

',1.6, 


0242  , 
.0825 


.2945 
02  32  , 


.  0780 
.3135 


0225  , 
.0765 


.  328  5 
0225, 


.0780 
.  3460 


0231 , 
.078  3 


-G 


,  .  2  5  , « 

1.0,1. 


.9,. 95 
1.7,1. 


.0254, 
, . 1005 


, .3367 
.0245, 


,  .  0960 
, .3566 


.02  42, 
, . 0967 


,.  3  7  70 
.0242, 


, .0946 
,  .  3  9  bu 


.0244, 
, .0953 


, .  40  75 
.  02  49, 


,  .09  73 
, . 4180 


3 ,.35, 
05,1.1 


,1.0,1 

8,1.9, 


.0284, 
, .1180 


,.3810 

.0273, 


,.1163 
,  .  4  lu5 


. 0268  , 
,.1156 


, . 4325 
.0265, 


. . 1150 
, .4550 


.0269, 

,.1152 


,  .  46  50 
.  02  74, 


,  .  1  1  6  3 
, .4800 


.2470 
.0300,, 
.C  735 
.  2565 


.0  394, . 
.  G  8  4  8  , 


.2690, 
. 6435,. 


.0910, 

.2825, 


2835 

299, 


0336,. 0389, .040 7, .04 38, 
. 0990, .  1135, .1333,  .1510 


•3060,. 3440, .4060, .4630 
0420, .0420,'.  0434, .0464, 


.  1050,.  12  13,  .  1425,  .  16  10 
•  3215..3  600 ,.4150, .47  00 


•  4  ,  .45 
,  1  •  15, 


.05,1. 

2.0,2. 


.0  3  20  , 
,.140  0 


,  .  4  4  9  5 

.0305  , 


,.  1  396 
,.4790 


.02  99  , 
,. 1 395 


,.  5040 
.0298, 


,.1395 
, .  52  8U 


.0300  , 
,  .  1400 


,  .  5590 
. 0  30 4  , 


.4860, .o 100 
0410, .0465  , 


,.5,. 55 

1  •  2  , 1 . 4 


1,1.2, 

2,2.4, 


.0360,. 
, • 1630, 


,  .  8  7  0  U  , 
.0347  ,  . 


,  .  16  4  0  , 
,.9000, 


.03  41  ,  . 

, . 1650, 


1 . 9200, 
.0  33  7,  . 


,  .  16  6  8, 

,.9350, 


•0336,. 
,  •  t  6  6  o  , 


,.  9400  , 
.0  343  ,. 


,  .  17  0  5, 
,.9900, 


38 
1  7 


,1.0500 

.0540,. 


0484, .0542, 
.1678, .1865 


•  52  00  » •  6 1>  5 0 
05  11  , .0580  , 


.  1 
.5 


♦.6, .65, 

,1.6, 


0420, 

.1850, 


1 .3100, 
0390  , 


.1873, 
1.  3400, 


0390, 

.  190  4, 


1.36  50, 
0  383  , 


.  1960 , 
1.4000, 


0386, 

.2000, 


1.4200, 
0396  , 


.2045, 

1.4700, 


1 8  , 

25, 


',  1.5200, 
062  3, 


0625, .0728, 
.2  130, .2400, 


1  .  15  00, 1.6200, 
0o7o*.0772  , 

25 

1  . 


Qi 


3 

5 


.  1062 
.  3075 


. 0496* 
.  1  130 
.3150 
.  0596, 


. 0500, 
.  I2b3 


.3170 
. 0586, 


.1300 
.  3320 


.  OS  7  3, 
13  30 


« . 1 25d 

« .  3  5  <  j  0 


.0583, 
,.  1328 
,  .3500 
.0586 


..1550 
,.  3  700 


•  J  5  8  3  , 
,.1555 
,.  3  700 
059 


,  •  55  C  u  » 
.0586,  . 


,.165  7 
.  5  J6 «' 


.  0  *,  9  7  »  • 
,  .  1795  , 
,.5210, 
0515  ,  . 


.1810, 
.51c  5, 


.  5  j-,0  , 
0  51  7 ,  * 


•  5o  00 
u  5d  0 » 


,  •  6  65  U 

5  , 


o  5  3  5  ,  . 
.  1995  , 
. 5/o J, 
■■  5  52  ,  . 


«  2  0 50  , 
.  5500  , 


0  tiuU  , 
•  2  1  0  0 
.  5  300 
06  2  2, 


.2260 
.  5950 


.  06‘*0  , 
,  .  £  52  0 
,.5850 


,.25  70 


,i.05  7 
.0765, 


,  •  2  6  2  6 
,.9100 


.0313, 
,  .  2  6  7  6 

,  •  5  0  00 
.  0  3  5  6  , 


»  »  2  7  1  u 
.  7  7o  C 


, . 3500 
.058 


,  .1505 
,  .  3550 


.0582, 
,  .  1575 


,  .  3600,  .395  3,  •  5  3  u  0  »  .56Go,.590u».65uO 
•<92, .052  5,  .(.'566,  .0653,  .0765,  .092C, 


,.1723, .I960, •  2  l  9  u  » .2500, .2650, .2860 
,.3800, .5100, .5600, . 5800,. 5150, .6700 


.0590, .0522,  .0575, .0668, .0  790, .0958, 
-  1800, .2050, .2230, .25 00, .2750, .3000 


•0903, 


, .2330, 
,1.3600, 


.  u  96  o  , 

, . 2900, 
,1.2600, 
. lulu  , 


, .2930, 

, 1.0 800, 


.  1 0  5  2  , 

, .2960, 


,  u  .  9  3  00  , 
.1098, 


, . 3090, 
,0.91500, 


. 1152  , 
,.3230, 


,.2550,. 3100, .3350, 
5500,0.6970,0.6750 


8  5 
29 


.70,24.95,26.  10,2  7. 
.  J  5 »  1 .14,  2 . b5,  4 . 


.55 ,12.  50 ,1  3.  45  ,  1  4  • 
.30,25.30,26.30,27. 


.92,  1.07,  3.10,  5. 
.  30, 17.50,  18.  75,20  . 


.20,28.20,28.20,28. 

4.93,  7. 


2o,2«.cO,24. 
17,  5.65,  7  . 


.  6  •  Gu  ,  l  7  • 
2  (  ,27.2>j,27. 


iv ,  7.56,10. 

J  -  ,  2  2  .  1 0,2  3  . 


2 0,28.2 9 ,28. 
65, l u . 6  7  ,13. 


,65,20.85 ,20.65,20. 
.25,  2-55,  5,30,  8, 


85 ,20.86,20 
10, 10.90,14 


20, 30. 
5  0,  8  • 


*,0,19. 

2  0,2  7. 


00,11. 

do , 25 . 


2  0 , 2d . 

7  0 , 1 6 . 


22. 

22. 


ds , 20. 
90, It. 


20, 31.10, 31.  10,  iH-l.o 
30,  4.10,  4.90,10.70, 


3 0, 2 0.  76, 22. J 0,2 3.20* 

2 0 , 2  7 . 2 0 , 2  7 . 20  ,  11*1.2 


25, 12. 5 0,13. dO, 15.  10, 
4u,2o.20,2b.20,2d.20, 


20, 2d. 20, 28. 20,  H=L.6 

00, Id. 20, 20. 45, 22.  70 


20.85,2 
2  O  .  8  5 , 2 


65 , 20. d 5, 20. 85 ,  11=2 .25 

90,14.90, 14.90, 14.90, 


CANARD  ORAG  POLAR 


l  NO  AO  4=11, 


I  A04X*  2  3, 
I  AG4Y  -12  , 


AT AB04--.22,-. 20,-. 18 
0.  ,  .02, . 04 • .06, 


.2, .4, .6, .8,1 


,  — .  16,  —  .14,“.  12,-. 10, -.08, -.06, -.04, -.02, 
.08,  .10,*  12, .14, .  16,  .18, .20, .22, 


.0,1. 2, 1.4, 1.6, 1.8, 2. 0,2. 2, 2.4, 


\  .1*0.2 


1  1  *  f  ■  .  4 


.0635,. 


.0011, 

.0016, 


.0262 


.0717, 

•0186, 


. OOlu , 
.0015, 


.0250, 


.0093,  .0057, 


00  32  , 


U  3  4<J  « 


.0057, 


.0428, 


0029, 


0324, 


\  1*0.6 


\  1*0.8 


l  3 
82 


.0010, 

.0014,  .0029, 


.0252,  .0362, 


553, 

159 


.0451, 


•  04  o  7  , 


.0512, 
.  0  0  8  6  , 


.0063, 


.05  12, 


. 0961  , 
.  009  V 


.0093, 


.  0  52  5 


.0  324, 
.0054, 


.0068  , 


.  U  5  3  3  , 


.0  362  , 
.0053, 


.  0086 


.0  719, 


.0616, 
.00  52  , 


.0  1  40  , 


35 


.0016  , 


.0196, 


.02  50, 
•  u029  , 


.01  32  , 


.0717  , 


.02  be  , 

.0029, 


.0015, 


.  0  i  8 


.1013, 


.0399, 

.0028, 


.0014, 


.0182, 


.0014 


.0009, 

.0014, 

,  .0028,  .0052,  .0094,  .0159,  .0255, 

.0399, 

1  1*1.0 

>  •  0  6  1 6  t  •  0  6  1  t  •  1  b  3  »  •  2  7  3  2  * 

t  1*1.2 


\  1=1.4 


1.181, 
.0317,  . 


.0014, 
.0017,  .0 


2  3,  . 06  76  , 


•  0  8  2  #  «  Oo'm  i 

.0255,  .0180, 


.OOld, 

.0022,  .0040, 

.0344,  •  04  *»  8  , 

•  1048,  .0  d6o  , 

. 0  J 1  5 ,  .0221  , 


.1552,  .3204, 


. o073, 
. 0565  , 

. 0702 , 

.0145, 


.0119, 
.0697 , 

.0555, 
.  0  0  8  7  , 


8 


1.18  1, 


.0344, 

.0040, 


.  o  1  8  0  , 
.0842  , 

•042b, 
•  0  4  t>  , 


.0022  , 


.0255, 

.  :  r.  I  4  , 


■  ■*  *  *  •  •  *  .  •  .  *  -  k  . 

■  J  •  s  j-  '  w  W.*  |'.  ■  - 
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1.212, 

2.656, 


1 . 2  85  , 


1 
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53 

95 
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1 
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1. 
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.993, 
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.936, 


1 .  i  a  2  , 
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Figure  24  —  Thrust  Required,  10,000  ft  (Case  1) 
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Figure  29  —  Sustained  Turn  Rate,  20,000  ft  (Cases  2&3) 
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Figure  30  —  Sustained  Turn  Rate,  30,000  ft  (Cases  2&3) 
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Figure  35  —  Specific  Excess  Power,  Sea  Level  (Cases  2&3) 
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